Nonequilibrium Stagnation-Line Radiative 
Heating for Fire II 


Christopher O. Johnston and Brian R. Hollis' 

NASA Langley Research Center, Hampton, VA, 23681 

Kenneth Sutton 1 

National Institute of Aerospace, Hampton, VA, 24060 


This paper presents a detailed analysis of the shock-layer radiative heating to the Fire II 
vehicle using a new air radiation model and a viscous shock-layer flowlield model. This new 
air radiation model contains the most up-to-date properties for modeling the atomic-line, 
atomic photoionization, molecular band, and non-Boltzmann processes. The applied viscous 
shock-layer flowlield analysis contains the same thermophysical properties and 
nonequilibrium models as the LAURA Navier-Stokes code. Radiation-flowfield coupling, or 
radiation cooling, is accounted for in detail in this study. It is shown to reduce the radiative 
heating by about 30% for the peak radiative heating points, while reducing the convective 
heating only slightly. A detailed review of past Fire II radiative heating studies is presented. 
It is observed that the scatter in the radiation predicted by these past studies is mostly a 
result of the different flowlield chemistry models and the treatment of the electronic state 
populations. The present predictions provide, on average throughout the trajectory, a better 
comparison with Fire II flight data than any previous study. The magnitude of the vacuum 
ultraviolet (VUV) contribution to the radiative flux is estimated from the calorimeter 
measurements. This is achieved using the radiometer measurements and the predicted 
convective heating. The VUV radiation predicted by the present model agrees well with the 
VUV contribution inferred from the Fire II calorimeter measurement, although only when 
radiation-flowfield coupling is accounted for. This agreement provides evidence that the 
present model accurately models the VUV radiation, which is shown to contribute 
significantly to the Fire II radiative heating. 
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Nomenclature 

Blackbody or Planck function written in terms of frequency (erg/cm 2 /s/sr/eV) 

velocity of light, equal to 2.997925xl0 10 cm/s 

electron charge, equal to 4.80298xl0" 10 cm 3/2 g 1/2 /s 

degeneracy for an atomic level i 

Planck’s constant, equal to 6.6256xl0' 27 erg-s 

wall-directed frequency-integrated radiative intensity (erg/s/cm 2 /sr) 

wall-directed frequency dependent radiative intensity (erg/cm 2 /s/sr/eV) 

frequency dependent emission coefficient (erg/cm 3 /sr) 

number density of an atom (particles/cm 3 ) 

ion number density (particles/cm 3 ) 

electron number density (particles/cm 3 ) 

wall-directed radiative heat flux, calculated in units of erg/s/cm 2 , but usually 
presented in units of W/cm 2 (lxlO' 7 W/cm 2 = lerg/s/cm 2 ) 
wall-directed frequency-dependent (in terms of eV) radiative heat flux 
convective heat flux (W/cm 2 ) 
time (s) 
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Tye 

T lr 

Z = 

aq r = 

A 

v = 

Subscripts 

EQ 

h v = 

Superscripts 

nc = 

sc = 


vibrational-electronic-electron temperature (K) 

translational-rotational temperature (K) 

distance along the stagnation line measured from the wall (cm) 

radiative flux absorbed by the beryllium calorimeter 

wavelength (cm) 

frequency (s' 1 ) 

assumes that Boltzmann and Saha equilibrium exists among the electronic levels 
indicates a spectral dependence in terms of eV 

indicates that the non-catalytic wall boundary condition was applied 
indicates that the super-catalytic wall boundary condition was applied 
indicates the radiative flux or intensity directed towards the vehicle wall 


Abbreviations 

eV =electron volts; the frequency in eV, labeled hv, is equal to 1.24xl0' 4 v '/ c 
VUV = vacuum ultraviolet; refers to the spectral region above 6 eV 


I. Introduction 

T he heat-flux to a vehicle at lunar-return conditions consists of a component due to the radiative flux from the 
high-temperature shock-layer gas 1 . The prediction of this radiative flux has been studied sporadically since the 
beginning of manned space-flight, with significant progress being made for NASA’s Apollo program. In support of 
this research, the Fire II flight experiment obtained radiative heating measurements at conditions similar to those for 
lunar-return 2 ' 3 . These measurements have since become a benchmark for new radiative heating prediction 
algorithms. The goal of this paper is to review the past predictions of these data and to present the results of a newly- 
developed air-radiation model. 

A review of past studies of the Fire II radiative heating is presented in Section II. The various predictions to the 
flight data are compared and the differences discussed. The flowfield and radiation model used by each researcher is 
defined and its influence on the comparison with the flight data is mentioned. The viscous-shock-layer flowfield 
model used in the present study is discussed in Section III. The number densities and temperatures from this model 
are compared with those from the LAURA Navier-Stokes code, and it is shown that these two flowfield models 
produce similar radiation predictions (when the same radiation code is applied to both). A brief description of the 
radiation model applied in the present study is presented in Section IV. The Fire II radiative and convective heating 
predicted using the present flowfield and radiation models are presented in Section V for the case where radiation- 
flowfield coupling is not considered. These results, which are known to be physically inaccurate, are presented to 
provide a baseline for comparing with the radiation-flowfield coupled cases. Section VI compares the radiation- 
flowfield coupled results of the present study with the flight data and various past predictions. The deviations of the 
present results from the flight data and previous studies are discussed, and the influence of non-Boltzmann radiation 
and flowfield coupling are assessed. The differences between the radiation predicted by a two-temperature chemical 
nonequilibrium VSL flowfield and a single temperature chemical equilibrium VSL flowfield are compared, which 
provides insight into the effect of multiple temperatures and chemical nonequilibrium on the radiative heating 
predictions and the effect of radiation-flowfield coupling. 

II. Review of Past Studies 

The Fire II flight experiment was flown in 1965 to obtain radiative heating data for a capsule with a scaled-down 
Apollo-shaped forebody 2 . Three non-ablating beryllium heat shields were attached to the vehicle, with the first two 
being jettisoned at selected points in the trajectory. Table 1 lists the flight conditions and wall temperature at various 
trajectory points along with the vehicle nose radius, which is different for the various heat-shields. Radiative heating 
data was obtained in the stagnation region with three different types of instrumentation. A total radiometer measured 
the radiative intensity behind a quartz window, which allowed the radiation in the spectral range of 0 to 6 eV to be 
measured. A spectral radiometer provided spectrally resolved intensity data in the spectral range of 2.2 to 4.1 eV 
with a resolution of about 0.04 eV and a reported uncertainty of ±23%. A beryllium calorimeter measured the 
convective heating plus the absorbed radiative flux. The absorbed radiative flux refers to the component that is 
absorbed by the beryllium calorimeter over the entire spectrum. The spectral absorptance of beryllium is reported by 
Cornette 3 . 
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Table 1. Trajectory points for Fire II cases 


t( s) 

Altitude 

(km) 

Density 

(kg/m 3 ) 

Velocity 

(km/s) 

Temp. 

(K) 

Wall 

Temp. (K) 

Physical 
R n (m) 

1634.0 

76.42 

3.72x1 0 s 

11.36 

195 

615 

0.935 

1636.0 

71.02 

8.57x1 0 -5 

11.31 

210 

810 

0.935 

1637.5 

67.05 

1.47 xIO -4 

11.25 

228 

1030 

0.935 

1640.5 

59.62 

3.86x1 0 -4 

10.97 

254 

1560 

0.935 

1643.0 

53.04 

7.80x1 0 -4 

10.48 

276 

640 

0.805 

1645.0 

48.37 

1 .32x10 3 

9.83 

285 

1520 

0.805 

1648.3 

41.60 

3.25x1 0 -3 

8.10 

267 

503 

0.702 


Many researchers have used the Fire II case as a benchmark for a new radiation code or flowfield 
solver 4,5 ' 6 ' 7,8 ' 9 ' 10 ' 11,12 ’ 13 . Table 2 lists those studies published since 1984, which are the most relevant to modern 
studies. Listed along with each researcher are the flowfield equations, chemistry model, radiation code, and method 
of obtaining the molecular and atomic state populations implemented in the particular study. All of the studies listed 
here include radiation-flowfield coupling, meaning the divergence of the radiative flux is included in the energy 
equation. 

Table 2. Summary of previous Fire II studies 


Researcher 

Ref. 

Flowfield 

Eqs/ 

Chemistry 

Radiation 

Code 

State 

Populations 70 

Sutton [1984] 

4 

Euler 

E 

RAD/EQUIL 

B 

Balakrishnan et al. [1985] 

5 

VSL 

E 

RAD/EQUIL 

B 

Gupta [1987] 

6 

VSL 

E 

RAD/EQUIL 

B 

Bird [1987] 

7 

DSMC 

NE 

custom 

NB 

Carslon [1989] 

8 

VSL 

E 

8-step model 

NB 

Park [1989] 

9 

VSL 

NE 

NEQAIR85 

NB 

Gaily [1991] 

10 

VSL 

NE 

RAD/EQUIL 

NB 

Greendyke et al. [1994] 

11 

NS 

NE 

LORAN 

NB 

Olynick et al. [1994] 

12 

NS 

NE 

NOVAR 

B 

Park [2004] 

13 

VSL 

NE 

NEQAIR04 

NB 


* 

VSL = Viscous Shock Layer, DSMC = Direct Simulation Monte Carlo, NS = Navier Stokes 

# E = Equilibrium, NE = Nonequilibrium 

% 

B = Boltzmann, NB = Non-Boltzmann 


The most widely analyzed data from the three measuring devices is the total radiometer data, which was 
presented in Figure 13 of Cauchon 2 as the radiative intensity integrated between 0 and 6 eV. Figure 1 compares the 
intensity values predicted by many of the researchers listed in Table 2 with the flight data. According to Cauchon 2 , 
the total radiometer data contains an uncertainty of + 20% , which is indicated by the error bars in the figure. The 
results of Olynick et al. 12 compare best with the data, while the results of Gupta 6 are also in good agreement 
throughout the trajectory. From Table 2, it is seen that Olynick et al. applied the code NOVAR, which assumed a 
Boltzmann distribution of electronic states, to a Navier-Stokes nonequilibrium flowfield. The NOVAR code is 
reported to be essentially the same as LORAN, except that it is configured for computational efficiency and assumes 
a Boltzmann distribution of electronic states, instead of applying the non-Boltzmann QSS model used in LORAN 
and NEQAIR. Note that Olynick et al.’s results compare well for the early trajectory points. This is unexpected 
since a large percentage of the shock layer is in nonequilibrium for these points, which should have caused the 
Boltzmann assumption to over-predict the radiation considerably 11 . In addition to the Boltzmann results, Olynick et 
al. also presented several cases that applied the non-Boltzmann QSS model. These results showed that the QSS and 
Boltzmann model agreed for the early nonequilibrium trajectory points and then diverged from each other for the 
later equilibrium points. This trend is inconsistent with other studies 14 , which show that a Boltzmann distribution of 
electronic states is approached in regions of chemical equilibrium. 

Another interesting feature of Figure 1 is the good agreement with the data obtained by Gupta 6 , who applied the 
RAD/EQUIL code to an equilibrium VSL flowfield. In a similar study, Sutton 4 also applied the RAD/EQUIL code 
to an equilibrium flowfield, although it was inviscid. Assuming that similar thermodynamic properties were used 
and that the radiation codes applied were actually the same, it woidd be expected that because of boundary layer 
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absorption, Gupta’s values should be slightly lower than Sutton’s. Surprisingly, Figure 1 shows the opposite trend 
throughout the trajectory. Furthermore, the results of Balakrishnan et al. 4 5 should agree very well with Gupta’s, since 
they both applied the RAD/EQUIL code to an equilibrium VSL flowfield (although their VSL techniques were 
different). It is seen in Figure 1 that Gupta’s and Balakrishnan et al.’s results do not compare well. Unlike Gupta’s 
values, Balakrishnan et al.’s values are lower than Sutton’s inviscid results, which indicate the expected influence of 
boundary layer absorption. 

The most state-of-the-art studies presented in Figure 1 are those of Greendyke et al." and Park 13 , which both 
applied a nonequilibrium radiation model to a coupled thermochemical nonequilibrium flowfield. Although these 
studies agree reasonably well with the data, they do not significantly improve upon the results of older studies. A 
reason for this lack of improvement, especially for the early trajectory points, is the uncertainty in the 
nonequilibrium flowfield modeling (which was not an issue with the older chemical-equilibrium studies). This point 
is made clear in a study by Hartung et al. 15 , who adjusted some of the unknown nonequilibrium flowfield parameters 
and studied their effect on the radiative heating for the Fire 1631, 1634, and 1637.5 second cases. It was shown that 
the integrated intensity, between 0 and 6 eV, varied for the 1634 and 1637.5 cases between 1.8 to 4.4 and 11.0 to 
20.0 W/cm 2 -sr, respectively, depending upon the nonequilibrium flowfield parameters chosen. Although 
recommended values for the parameters are provided in the literature 16 ’ 17 ’ 18 , the evidence supporting these values is 
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Figure 1. Comparison of past predictions for the partial intensity at the wall with the Fire II total radiometer data 

throughout the trajectory. 

The frequency-integrated data from the spectral radiometer (between 2.2 and 4.1 eV), presented in Figure 13 of 
Cauchon 2 , has been the subject of comparison by only a few researchers. The comparison between the data and 
predictions are shown in Figure 2. The two “Flight Data” lines represent the upper and lower limit of the data 
scatter. The predictions all agree reasonably well, although none of them stay within the data scatter throughout the 
entire trajectory. An interesting aspect of this comparison is the close agreement of the values predicted by Sutton 
and Balakrishnan et al., as opposed to Figure 1 where the two differ significantly, especially near the peak heating 
point. A few studies also compared their predictions with the spectrally-resolved data 19 from the spectral 
radiometer 1 ’ 15 . Park 13 showed that his predicted spectrum compared well with the measured spectrum for the 1643 s 
point. This result is consistent with the decent comparison of the frequency-integrated values presented in Figure 2. 
Hartung et al. 15 compared their computed spectra obtained with the LORAN code, with the data for the 1631, 1634, 
and 1637.5 s trajectory points. The agreement with the data for these nonequilibrium cases was not very good, but 
since they did not present values for the integrated intensity between 2.2 and 4.1 eV, it is difficult to quantify the 
disagreement. From Figure 5 of Hartung 20 and Figures 6 and 7 presented by Greendyke 21 , it is suspected that the 
non-Boltzmann modeling of N 2 and N 2 + are likely responsible for this poor agreement with the data. The 2.2 to 4.1 
eV region of the spectrum is dominated by the N 2 (l + ) and N 2 + (H) bands, which emit strongly in the nonequilibrium 
region near the shock. The LORAN code applied by Hartung uses the same non-Boltzmann QSS model as 
NEQAIR. 

The Fire II calorimeter data presented by Cornette’ provide heating values that contain the convective heating 
plus the contribution of the radiative flux absorbed by the beryllium calorimeter. A comparison of the calorimeter 
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data with various predictions, throughout the trajectory, is presented in Figure 3. The +5% error bars on the data in 
this figure indicate the measurement uncertainty according to Cornette’. Although the flight data cannot be separated 
into convective and radiative components, Figure 4 compares the calculated values of these components, whose 
values add together to produce the values shown in Figure 3. The best comparisons with the data in Figure 3 are the 
results of Olynick et al. 12 , as was the case for the radiometer data in Figure 1. It is also seen in Figure 3 that the 
results of Sutton 4 and Gupta 5 6 are in close agreement with each other throughout the trajectory. This agreement is due 
to their offsetting differences in the radiative and convective components. Surprisingly, Sutton’s radiative 
component is larger than Gupta’s, which is the opposite of what was seen in Figure 1 for the 0 to 6 eV integrated 
intensity. This indicates that Sutton predicted a larger radiative contribution from the spectral region above 6 eV 
(the vacuum ultraviolet) than Gupta. This difference is likely due to boundary layer absorption, which influences the 
vacuum ultraviolet more than any other spectral region. The results of Greendyke et al." provide the worst 
comparison with the data. It is seen that both the convective and radiative predictions are too low, especially near 
peak heating. Note that the large difference in convective heating predicted by Greendyke et al. and Olynick et al. is 
present even though both studies applied a nonequilibrium Navier-Stokes flowfield with a super-catalytic wall 
boundary condition. 



Figure 2. Comparison of past predictions with the frequency-integrated 
Fire II spectral radiometer data. 



Figure 3. Comparison of past predictions for the total absorbed heat ilux 
with the Fire II calorimeter data. 
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Figure 4. Comparisons of the predictions by various researchers for the radiation absorbed 
by the calorimeter (aq,') and the convective heating (q,). 

III. Flowfield Model 

The flowfield and radiation model applied in the present study is discussed in detail by Johnston 22 . A viscous 
shock layer (VSL) method 2 ' was applied for the flowfield model. This method allows for a computationally efficient 
solution that retains all of the important flowfield physics for the stagnation region of a blunt body. Only the 
stagnation line flowfield was considered in the present study, although the same method may be used to obtain 
solutions downstream of the stagnation line. The same kinetic models and thermophysical properties are applied in 
the present VSL code as are applied in the LAURA code 24 . The air model applied consists of the following 11- 
species: N 2 , N 2 + , 0 2 , 0 2 + , NO, NO + , N, N + , O, 0 + , and e". The forward chemical rates were taken from Park 25 , while 
the backward rates were obtained by applying detailed balancing; with the equilibrium constants calculated using 
curve-fits from McBride et al. 26 . These curve-fits were also used to calculate the specific heat and enthalpy values 
for each species. The collision cross-sections required for the calculation of the diffusion, viscosity, and thermal 
conductivity coefficients were taken from Gupta et al. 27 . Multicomponent diffusion was modeled using the 
“approximate-corrected” approach presented by Sutton and Gnoffo 28 . 

To validate the developed VSL model for lunar return conditions, comparisons were made with Navier-Stokes 
results produced by the LAURA code for Fire II cases. The temperature and species number densities (relevant to 
the radiation calculation) predicted by the VSL method and the LAURA code are compared in Figure 5 for the Fire 
II 1636 s trajectory point. This point was chosen because the shock-layer contains both a region of chemical 
equilibrium and nonequilibrium. In these figures, the vehicle wall is located at z = 0, where z is the distance along 
the stagnation line. The consequence of the different treatments of the bow shock wave by the two methods is 
clearly shown in these figures. The “discrete shock” model of the VSL approach is seen to result in a larger 
translational-rotational temperature directly behind the shock, even with the application of the shock-slip equations. 
Along with beginning at a higher temperature directly behind the shock, the T tr values in Figure 5 are seen to behave 
slightly different in the nonequilibrium region directly behind the shock. This is another consequence of the discrete 
shock treatment, which results from T tr beginning its relaxation process behind the shock at a higher temperature 
than for the LAURA case. The T ve values, on the other hand, are in relatively good agreement in the nonequilibrium 
regions. This is fortunate because T ve governs the radiation. Other than the differences near the shock, the 
temperature profiles shown in Figure 5 closely agree throughout the rest of the shock-layer. The locations at which 
the two temperatures equilibrate compare well between the methods, as do the values of the temperatures in the 
equilibrium regions of the layer, which agree within 1% for this case. This close agreement in the equilibrium region 
of the shock-layer is required for the prediction of similar radiative emission values, because of the exponential T ve 
dependence of the radiative emission. This exponential dependence is due to the Boltzmann distribution of the upper 
electronic state of the various radiative transitions, and hence is only true in the regions of chemical equilibrium 
where a Boltzmann distribution is approached. Conversely to this exponential dependence on temperature, most of 
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the radiative emission depends linearly on the species number densities (the atomic bound-free emission has 
essentially a quadratic dependence). The disagreement in the equilibrium regions of no more than 5% for the present 
model is therefore regarded as sufficient. 



Figure 5. Comparison of the VSL and LAURA stagnation line temperature profiles and 
number densities for the Fire 1636 case. 

To make certain that the VSL method provides a sufficient flowfield model for predicting radiative heating, the 
radiation model to be discussed in the next section was applied to the Fire II 1634 and 1636 s flowfields. The 
radiation model for these calculations includes the non-Boltzmann model for the atomic and molecular electronic 
states. Figure 6 presents the wall-directed radiative flux for these cases resulting from the uncoupled VSL and 
LAURA flowfields. The values at the wall (z = 0) agree within 5% for the 1636 s case and 1% for the 1634 s case, 
while the values throughout the shock-layer agree closely as well. This agreement throughout the shock layer is 
required for the accurate modeling of radiation-flowfield coupling, which although not included here, is dependent 
on the wall-directed (and shock-directed) radiative flux at each point through the shock-layer, as will be discussed in 
Section VI. 



Figure 6. Wall directed radiative heat flux values predicted by the non-Boltzmann radiation 
models for the various flowfields. 
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IV. Radiation Model 

The radiation model applied in the present study is discussed in detail by Johnston" 2 . This model is based on a set 
of atomic levels and lines obtained from the National Institute of Standards and Technology (NIST) online 
database 29 and atomic bound-free cross sections from the Opacity project’s online TOPbase 30 The molecular band 
systems are treated using a smeared-rotational band (SRB) model’ 1 . The molecular data for modeling these band 
systems are obtained from Laux 32 , except for the vacuum ultraviolet N 2 systems, which are obtained from various 
other sources 33,34,35 . The non-Boltzmann modeling of the atomic and molecular electronic states is based on a set of 
electron-impact excitation rates compiled from the literature and presented in detail by Johnston 22 . Following the 
work of Park 36 , the quasi-steady state assumption is made when solving the Master Equation. The tangent-slab 
approximation is applied to calculate radiative flux and the divergence of the radiative flux, which is required for the 
radiation-flowfield coupling procedure. 

V. Uncoupled Radiative Heating for Fire II 

The Fire II flight experiment was discussed in detail in Section II. Along with presenting the flight data, it was 
shown in this section that there is a large scatter among the previous theoretical predictions of this data. 
Furthermore, no prediction method compared well throughout both the early nonequilibrium region of the trajectory 
and in the later mostly equilibrium peak-heating region. It was shown that the more recent and sophisticated 
prediction methods did not lead to a consistently better comparison with the data relative to older prediction 
methods. The goal of this section and Section VI is to compare the results of the new flowfield and radiation models 
presented in Sections III and IV with the flight data and previous predictions. It will be shown that the present model 
provides a better comparison with the data, on average over the trajectory, than any previous study. Although the 
uncoupled results are known to be physically inaccurate, they are presented in this section so that they may be 
compared with the uncoupled values predicted by previous studies. This allows the differences in the flowfield and 
radiation models to be examined without the complication of radiation-flowfield coupling. Furthermore, it will be 
insightful to compare these uncoupled values with the coupled values presented in the next section. 

The chemical kinetics and thermophysical properties applied in the present flowfield model were discussed in 
Section III. Examples of uncoupled Fire II flowfields were also presented in that section, and shown to provide 
uncoupled radiative heating values consistent with those resulting from a LAURA (Navier-Stokes) flowfield. Table 
3 lists the various radiative components predicted with the present uncoupled nonequilibrium VSL model and 
radiation model discussed in Sections III and IV. The q and I w values listed here represent the wall-directed 
radiative flux and wall-directed intensity, respectively, at the wall (z = 0). The numbers in parenthesis in the header 
column represent the spectral range, in eV, included in each value. The aq,~ values represent the radiation absorbed 
by the beryllium calorimeter, which is a quantity required for comparison with the calorimeter data (as mentioned 
previously in the discussion of Figure 4). Also required for the comparison with the calorimeter data is the 
convective heating ( q c ), which is listed in this table assuming both a super-catalytic (sc) and non-catalytic (nc) wall. 
The difference in the radiation predicted by these two catalyticity assumptions is negligible, and so only the results 
from the super-catalytic case are presented here. 


Table 3. Uncoupled convective and radiative heat flux values (W/cm 2 ) and radiative intensities (W/cm 2 -sr) 
for the Fire II vehicle obtained using the NE VSL flowfield model. 


t (s) 

sc 

nc 

ac/r 

( 0 - 18 ) 

ccq r ' 
(0-6) 

aq,~ 

(6-18) 

9r 

( 0 - 18 ) 

dr 

(0-6) 

dr 

(6-18) 

h 

( 0 - 18 ) 

L 

(0-6) 

I w 

(6-18) 

1634.0 

221 

115 

19.2 

5.53 

13.7 

25.7 

10.9 

14.8 

4.9 

1.8 

3.2 

1636.0 

316 

195 

71.6 

18.7 

52.9 

96.1 

38.5 

57.7 

19.5 

6.4 

13.0 

1637.5 

398 

266 

152 

40.8 

111 

206.3 

85.1 

121.2 

42.6 

14.5 

28.1 

1640.5 

578 

439 

427 

140 

287 

615.5 

297.0 

318.5 

129.6 

51.8 

77.8 

1643.0 

777 

689 

500 

197 

303 

755.5 

415.8 

339.7 

156.2 

71.9 

84.4 

1645.0 

835 

743 

253 

124 

129 

403.7 

256.5 

147.2 

79.1 

42.9 

36.3 

1648.3 

756 

670 

25.6 

13.2 

12.4 

41.8 

25.2 

16.6 

7.1 

4.0 

3.1 


The two previous studies that present uncoupled results over the majority of the trajectory are the studies by 
Olynick et al. 12 and Greendyke et al. 11 . The comparison between the q,~ and/„(0-6) values presented in Table 3 and 
those presented in Table V of Olynick et al. indicates agreement within 10% over most of the trajectory. The values 
presented by Greendyke et al. in their Figure 4 and 5, on the other hand, are lower than the present values by up to 
30%. The fact that Olynick et al.’s values agree well, while Greendyke et al.’s do not, is interesting considering that 
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Olynick et al. applied a Boltzmann model and Greendyke et al. applied a non-Boltzmann model for the electronic 
state populations. As mentioned previously, the present study also applies a non-Boltzmann model, which is known 
to result in less radiation than the Boltzmann model, especially in regions of significant chemical nonequilibrium. 
Note that the flowfields predicted by Greendyke et al. and Olynick et al. were significantly different (see Figure 3 of 
Olynick et al. 12 ), even though both were nonequilibrium Navier-Stokes solutions. 

After the completion of the present research, the study by Hash et al. 37 of the uncoupled convective and radiative 
heating for Fire II became available to the authors. This work compared the convective heating predicted by various 
state-of-the-art Navier-Stokes codes. Because the present flowfield model applies the same energy equation and 
thermophysical modeling as the LAURA code, the LAURA results presented by Hash et al. should compare well 
with the present results. The super-catalytic and non-catalytic convective heating predicted by the LAURA code at 
1636 s agree with the present results within 2%. For the 1643 s case, the super-catalytic value agrees within 3% 
while the non-catalytic value agrees within 9%. The current radiation results are 26% greater, in the 0-6 eV range, 
and 47% greater, in the 6-18 eV range, than the NEQAIR results presented by Hash et al. for the 1636 s case. 

VI. Coupled Radiative Heating for Fire II 

The radiation components predicted by the present radiation-coupled flowfield model are presented in Table 4 
for the Fire II case. The column definitions are identical to those of Table 3. To illustrate the differences between 
these values and the uncoupled results discussed previously, Figure 7 compares q c , q r ~, and aq,' between the coupled 
and uncoupled cases. It is seen that the q,~ and aq,' values are reduced by 20 - 30% throughout the main heating 
portion of the trajectory. This reduction is slightly larger than the 15 - 25 % reduction reported by Olynick et al. 12 
and the 5 - 18% reduction reported by Greendyke et al. 11 . The smaller reduction found by Greendyke et al. is most 
likely a result of their significantly smaller uncoupled radiation prediction, which reduces the coupling effect. The 
discrepancy with Olynick et al.’s result is likely a result of differences in the flowfield modeling. A significant 
difference between Olynick et al.’s flowfield model and the present model is that Olynick et al. obtains the 
vibrational-electronic-electron temperature by solving the vibrational energy equation, and not the vibrational- 
electronic-electron energy equation, as is done in the present model. The consequence of this treatment by Olynick 
et al. is that the divergence of the radiative flux does not appear in the vibrational energy equation, even though the 
temperature obtained from this equation governs the magnitude of the radiation. The influence of radiation on the 
behavior of the two temperatures is therefore fundamentally different for this model because radiation does not 
directly influence the vibrational-electronic-electron temperature. For the high-temperature shock-layers of present 
interest, the vibrational energy is very small throughout most of the shock-layer because the majority of the 
molecules are dissociated. Furthermore, the coupling of the vibrational energy mode to the translational energy 
mode is represented through vibrational-translational energy relaxation, which is shown in Eq. (55) of Gnoffo et al. 24 
to be proportional to the mass fractions of the molecules. For the present highly-dissociated conditions, this term is 
small, which results in the vibrational energy mode being weakly coupled to the translational energy mode. This 
indicates the possible inadequacy in solving the vibrational energy equation instead of the vibrational-electronic- 
electron energy equation for the high-temperature shock-layers of present interest. 


Table 4. Coupled convective and radiative heat flux values (W/cm 2 ) and radiative intensities (W/em 2 -sr) for the Fire II vehicle 
obtained using the NE VSL flowfield model. 


t (s) 

q c 

sc 

<lc 

nc 

cay," 

(0-18) 

aq r ~ 

(0-6) 

aq,: 

(6-18) 

<lr 

(0-18) 

<lr 

(0-6) 

q r 

(6-18) 

4 

(0-18) 

4 

(0-6) 

4 

(6-18) 

1634.0 

217 

113 

17.2 

5.3 

11.9 

23.2 

10.3 

12.9 

4.4 

1.7 

2.8 

1636.0 

304 

187 

59.4 

16.7 

42.7 

80.7 

34.2 

46.5 

16.2 

5.7 

10.5 

1637.5 

377 

253 

119 

35.0 

83.9 

165 

72.6 

92.0 

33.6 

12.3 

21.3 

1640.5 

535 

410 

309 

112 

196 

455 

236 

219 

94.2 

40.8 

53.4 

1643.0 

733 

654 

370 

159 

211 

567 

334 

233 

115 

57.3 

57.8 

1645.0 

807 

721 

208 

107 

101 

338 

223 

115 

65.2 

37.0 

28.3 

1648.3 

752 

672 

25.6 

13.2 

12.4 

41.7 

25.1 

16.6 

7.1 

4.0 

3.1 


The influence of radiation-flowfield coupling on a two-temperature chemical nonequilibrium flowfield should be 
similar to that predicted by a single-temperature chemical equilibrium model, assuming that the electron- 
translational energy exchange term is properly treated (for the two-temperature model). By comparing the present 
two-temperature chemical nonequilibrium VSL flowfield (NE VSL) with a chemical equilibrium VSL flowfield (E 
VSL), the influence of coupling on the two models may be compared. Figure 8(a) compares the partial intensity (0 - 
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6 eV) at the wall predicted by the uncoupled and coupled NE VSL and E VSL predictions. The decrease in intensity 
with the addition of coupling is seen to be very similar for both flowfield models. This confirms the previous 
statement regarding the similarity of the coupling effect for the two flowfield models. Figure 8(b) presents the 
coupled and uncoupled temperature profiles for the NE VSL and E VSL models at the 1643 s condition. Note that 
the vertical-axis of this figure ranges from 10,000 K to 12,000 K, for clarity. Except for the differences in the 
boundary layer, which is a result of the different diffusion modeling and possible nonequilibrium effects, the NE 
VSL and E VSL models compare very well for both the coupled and uncoupled cases. The slight separation seen 
between the models throughout the layer is only 10 to 20 K, which is negligible relative to the temperatures of 
roughly 11,000 K. 



Figure 7. Uncoupled and coupled values for the convective and radiative heat flux along the Fire II trajectory. 
The aq, values represent the radiative flux absorbed bv the calorimeter. 

4 

x 10 



t(s) 



Figure 8. a) Partial wall intensity values (0-6 eV) predicted by the coupled and uncoupled application of the NE and E VSL flowfield 
models, b) Uncoupled and coupled stagnation-line temperature profiles for t = 1643 s. 


To assess the influence of non-Boltzmann modeling and radiation- flowfield coupling on the radiative flux 
throughout the shock layer, and on the flux that ultimately reaches the wall. Figures 9 and 10 present the wall- 
directed spectrally-integrated radiative flux profiles through the shock-layer and the radiative flux spectrum at the 
wall. In these figures, both the non-Boltzmann and Boltzmann results are presented for the coupled case, while only 
the non-Boltzmann result is presented for the uncoupled case. The non-Boltzmann and Boltzmann terminology used 
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4 r (W/cm ) q- (\y/ cm 2 ) 


here specifies that the atomic and molecular electronic states are modeled by their respective Boltzmann or non- 
Boltzmann models. These models are discussed in detail in Johnston’ 2 . The accumulated flux values presented in 
Figures 9 and 10 allow for the contributions from the various radiative mechanisms to be compared (note that the 
spectmm for only the coupled non-Boltzmann case is presented in each figure). For the 1636 s case, the non- 
Boltzmann influence is seen to be large. This is a result of the relatively large region of nonequilibrium directly 
behind the shock, which was shown previously in Figure 5. From Figure 9, the increased flux for the Boltzmann 
case is seen to come from the atomic lines between 1.3 and 1.8 eV and the N 2 + first-negative band system located 
between 2.0 and 5.0 eV. The atomic bound-free continuum located above 12.0 eV also provides some excess 
radiation for the Boltzmann case. For the 1643 s case, Figure 10 shows that the Boltzmann influence is small. This is 
a result of most of the shock-layer being in equilibrium, which was shown by the temperature profiles presented in 
Figure 8(b). The influence of coupling for both the 1636 and 1643 s cases is seen to result mostly from the atomic 
lines located between 1.3 and 1.8 eV and the atomic bound-free continuum located above 12.0 eV. 



Figure 9. Wall-directed radiative flux profiles and wall spectra resulting from the coupled and uncoupled cases for t = 1636 s. 



Figure 10. Wall-directed radiative flux profiles and wall spectra resulting from the coupled and uncoupled cases for t = 1643 s. 


The Fire II total radiometer data and the theoretical predictions of this data by past studies were presented 
previously in Section II. This data is compared in Figure 11 with the present radiation model coupled to the 
nonequilibrium (NE VSL) and equilibrium flowfield (E VSL) models. The predictions from past studies are also 
shown in these figures for comparison. Figure 11(b) shows the values for all trajectory points considered, while 
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Figure 11(a) focuses on the three earliest points, which are important because they contain a significant amount of 
nonequilibrium radiation. From Figure 11(b) the differences between the flight data, previous predictions, and the 
present predictions can clearly be seen for the 1643 s point. For this trajectory point, the study by Olynick et al. 12 is 
the only previous prediction that is closer to the data than the present NE VSL result. As discussed previously in this 
section, the Olynick et al. values are likely larger because of the reduced influence of radiation-flowfield coupling, 
resulting from the different treatment of the energy equations and temperatures. The Boltzmann assumption applied 
by Olynick et al. should not influence the radiation at 1643 s because the flowfield is mostly in equilibrium, as 
indicated in Figure 10. Note that the present equilibrium values (E VSL) are slightly larger than the values predicted 
by Gupta 6 , who also applied an equilibrium VSL method. Assuming that the flowfield models were nearly 
equivalent, this trend is consistent with the fact that the present radiation model predicts larger radiation values for a 
given condition than the RAD/EQUIL code applied by Gupta. Also, the equivalent nose radius applied by Gupta 
(and Sutton) is 18% smaller than the currently applied radius, which would have possibly reduced the shock- 
standoff distances, and therefore the radiation values, below those predicted by the present model. Gupta does not 
present his predicted shock-standoff values, so this cannot be confirmed. The reasonable agreement with Park’s 13 
prediction at 1643 s is noted. It should also be noted that the under-prediction at the 1637.5 point seen in Figure 
11(a) is possibly a result of the radiometer window reaching its melting point, which would have caused the data at 
this trajectory point to be inaccurate. Of all the flight data values shown in Figure 11, this is the only point that was 
not considered a “prime” data period by Cauchon 2 . For the early nonequilibrium points presented in Figure 11(a), 
the present results provide the most consistent agreement with the data, excluding the results of Olynick et al. 12 . The 
assumption made by Olynick et al. of a Boltzmann distribution of the electronic state populations should have 
caused a significant over-prediction of the radiation at the early trajectory points, as indicated by the behavior shown 
in Figure 9. Since this behavior was not predicted by Olynick et al., it may be concluded that their flowfield model 
was inconsistent with the present nonequilibrium VSL method, as well as past and present results from the LAURA 
Navier-Stokes code. A surprising result of Figure 1 1(a) is the relatively close agreement between the NE VSL and E 
VSL models, even though the shock layer is in considerable chemical nonequilibrium. This is a result of the non- 
Boltzmann model suppressing the radiation from the regions of chemical nonequilibrium in the NE VSL flowfield, 
as shown in Figure 9. Note that if the comparison between the NE and E VSL flowfields was made using a 
Boltzmann radiation model, the NE VSL prediction would be significantly larger than the E VSL prediction. Flence, 
if nonequilibrium chemistry is accounted for, a non-Boltzmann model must be applied. 



t (s) t (s) 


Figure 11. Comparison of the measured partial intensity (0-6 eV) with various predictions, including the 
results of the present radiation model applied to the NE and E VSL codes: a) Early nonequilibrium trajectory points; 

b) Peak heating trajectory points. 


The Fire II spectral radiometer data, which measured the intensity in the 2.2 - 4.1 eV spectral range, was 
presented and compared with past predictions in Section II. The results of the present radiation model, coupled to 
the nonequilibrium (NE VSL) and equilibrium flowfield (E VSL) models, are compared with this data in Figure 12. 
The flight data is represented by two lines, representing the upper and lower limit of the data, which contained 
significant scatter. Good agreement is seen for the NE VSL method throughout the entire trajectory, while the E 
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VSL results are slightly lower than the data. The radiation contained in this limited spectral range is due mostly to 
the N 2 4 first negative band system. For the 1634 to 1640 s points, this radiation is controlled significantly by the 
rates chosen by Johnston 2 " for the excitation of the N 2 + (B) state. The agreement of the present model with the data at 
these points provides confidence that appropriate rates were chosen. Note that both Sutton’s and Gupta’s values are 
larger than the present predictions for the intensity from this 2.2 - 4.1 eV spectral range. Since Figure 1 1(b) shows 
that the contribution from the 0 - 6 eV range in these studies is actually smaller than the present predictions, then the 
radiation from the 0 - 2.2 eV range must have been significantly under-predicted by Sutton and Gupta relative to the 
present model. This is likely a result of the increased line radiation in the 1 - 2 eV range resulting from the updated 
set of the atomic lines applied in the present model, and discussed by Johnston" 2 . 



1634 1636 1638 1640 1642 1644 1646 1648 


t(s) 

Figure 12. Comparison of the partial intensity (2.2 - 4.1 eV) measured by the Fire II spectral 
radiometer with the present coupled results and previous predictions. 

The Fire II total calorimeter measurements presented by Cornette 3 contain the convective heating component 
along with the radiation absorbed by the beryllium calorimeter. This measured quantity, which will be labeled q tota i, 
may be written as 

<hotal=<lc +a <lr~ (!) 

where q c is the convective heating and aq r represents the radiation absorbed by the beryllium calorimeter. The 
absorbed radiation was obtained by multiplying the frequency-dependent radiative flux by the frequency-dependent 
absorptance (of polished beryllium) presented by Cauchon". The two components of Eq. (1) are presented in Tables 
3 and 4, with both the super-catalytic (sc) and non-catalytic (nc) values presented in each table. Figure 13 compares 
the super-catalytic and non-catalytic values for q, ota / with the flight data. For both the coupled and uncoupled super- 
catalytic cases, considering the flux from the 0-18 eV, a consistent over-prediction of flight data is seen to exist at 
all trajectory points except the last two. Interestingly, the coupled non-catalytic case, considering the flux from the 0 
-18 eV, is seen to compare well for all trajectory points, except the 1645 s point. This comparison supports the 
argument that the beryllium surface of the calorimeter, at temperatures ranging from 600 to 1,600 K, is closer to 
being non-catalytic than super-catalytic. Note that, as indicated by Figure 3, the poor comparison at 1645 s is 
consistent with all previous theoretical predictions, and has not been explained. Although, due to the radiometer 
window reaching its melting point at this trajectory point, the radiometer measurement is not accurate enough to 
confirm the radiation prediction between 0 and 6 eV. Therefore, it is possible that some unknown radiation 
mechanism, not predicted by any previous study, contributed to the apparent shift of the peak heating point. Note 
that Sutton 4 showed that this shift was not due to the values used for the beryllium absorptance. 

Due to the strong absorption of VUV radiation and the limitations of instrumentation windows to transmit in this 
spectral region, it is difficult to obtain quality data to assess the accuracy of theoretical predictions for VUV 
radiation. The Fire II calorimeter provides the best flight experimental data to infer the magnitude of VUV radiation, 
although it requires an assumption that the predictions for the convective heating and the radiation in the 0 to 6 eV 
spectral range are reasonably accurate. To indicate the influence of the vacuum ultraviolet (VUV) radiation in the 
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present calculations, Figure 13 presents the super-catalytic convective heating plus the absorbed radiative heating 
between 0 and 6 eV for both the coupled and uncoupled cases. The values presented here were taken from Tables 3 
and 4, which means that the coupled cases considered the entire 0 to 18 eV radiation contribution while calculating 
the flowfield (even for the cases where only the 0-6 eV contribution is assumed absorbed by the calorimeter). It is 
seen that both the uncoupled and coupled values are below the flight data at the peak heating points. 



Figure 13. Comparison of the present results, assuming a catalytic and non-catalytic wall, with the Fire II calorimeter data. These are 

coupled values from the nonequilibrium VSL flowfield. 

While Figure 13 implied the importance of the VUV radiation component by showing that it was required to 
achieve a good comparison between the calorimeter data and predictions, a stronger conclusion may be made 
regarding the VUV radiation if the comparison between the radiometer data and predictions are included in the 
discussion. Table 5 achieves this by converting the intensity (between 0 and 6 eV) measured by radiometer 
(Av,mms(0-6)) to the radiative flux absorbed by the calorimeter between 0 and 6 eV (aq~ r j n f erre j(0-6)), which is then 
subtracted along with the predicted convective heating from the actual calorimeter measurement q t „t a i, meas , resulting 
in the radiative flux absorbed by the calorimeter between 6 and 18 eV (aq r ,mferr<d£>- 18)). The subscript “inferred” is 
used to indicate that the value is a result of the actual calorimeter and radiometer data, along with the predicted 
convective heating and the ratio between the predicted aq r ~(0-6) and /„.( 0-6) values, labeled /3 in the table. Note that 
the predicted magnitude of the VUV radiation enters this analysis only indirectly through its influence on the 
coupled-radiation flowfield used to calculate q c , aq r ~i„f erred { 0-6), and 1J0-6). This makes the comparison between the 
“inferred” and predicted values in the last two columns of Table 5 meaningful. It is seen that, except for the 1643 
and 1645 s points, the predicted values of the absorbed VUV radiation (a<7,'(6-18)) fall between the ca/).;„/ e ,. re ^6-18) 
values obtained using the super-catalytic (nor in parenthesis) and non-catalytic (in parenthesis) q c values. The 
negative aq~ r , inferred 0-18) and aq~ r , inferred 6-18) values in this table are a result of the q to mi,meas being lower than the 
predicted super-catalytic convective heating. The non-physical nature of these negative values suggest that the heat- 
shield for this case is more accurately modeled with the non-catalytic assumption, which are shown in Table 5 to 
result in positive aq' rin f erre j( 0-18) and aq~ rin f e „. ed { 6-18) values throughout the trajectory. The disagreement at 1645 s 
between the predicted and inferred values is not surprising, since it was not considered a “prime” data point by 
Cauchon 2 (although 1637.5 and 1640.5 were not either). On the other hand, the over-prediction of aq' r ( 6-18) relative 
to aq ri ,f elTe J^ 6-18) (both the super-catalytic and non-catalytic value) at 1643 s is concerning because it is both a 
prime data point and the peak radiative heating point. A possible reason for this over-prediction is that the tangent- 
slab approximation, applied for the radiation transport calculation, does not account for the shock-layer curvature. 
Past studies concerning the accuracy of the tangent slab approximation have concluded that a detailed three- 
dimensional transport calculation produces a flux value about 10 - 20% lower than that predicted by the tangent slab 
approximation' 8 . Table 6 presents the values defined in Table 5, except now with both aq r (Q-6) and aq',( 6-18) 
multiplied by a factor of 0.85. Only the values that are influenced by this factor are presented in Table 6, with aq 
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,(0-6), aq',{0-6), and aq rtin f erre f 0-6) each being reduced by 0.85 while aq r j n f erre f6-18) is increased due to the 
reduction in aq' rM f erre f 0-6). With this adjustment, the 1643 s prediction for aq',{ 6-18) now falls between the super- 
catalytic and non-catalytic “inferred” values. Note that the intensity’ predictions do not require the tangent-slab 
approximation, and so the previous conclusions regarding the radiometer data are not influenced by this discussion. 


Table 5. Manipulation of the Fire II radiometer and calorimeter data with predictions to isolate 
the absorbed flux measured by the calorimeter in the 6 - 18 eV spectral range. 

The columns are defined below the table. 


f(s) 

1 w, meas 

(0-6) 

P 

inferred 

(0-6) 

ary," 

(0-6) 

Q total, meas 

inferred 

(0-18) 

aq,: 

(0-18) 

inferred 

(6-18) 

aq,: 

(6-18) 

1634.0 

1.3 

3.1 

4.1 

5.3 

175.0 

-42.0 (62.0) 

17.2 

-46.1 (57.9) 

11.9 

1636.0 

5.2 

2.9 

15.2 

16.7 

290.0 

-14.0 (103) 

59.4 

-29.2 (87.8) 

42.7 

1637.5 

13.0 

2.8 

37.0 

35.0 

410.0 

33.0 (157) 

119 

-4.0 (120) 

83.9 

1640.5 

35.0 

2.7 

96.1 

112 

750.0 

215 (340) 

309 

119 (244) 

196 

1643.0 

64.0 

2.8 

177.6 

159 

1025.0 

292 (371) 

370 

114 (193) 

211 

1645.0 

26.0 

2.9 

75.2 

107 

1130.0 

323 (409) 

208 

248 (334) 

101 

1648.3 

5.0 

3.3 

16.5 

13.2 

750.0 

-2.0 (78.0) 

25.6 

-18.5 (61.5) 

12.4 


Iw.meas (0-6) = Intensity measured by the radiometer (W/cm 2 /sr) 

P= «g,’(0-6)// vl ,(0-6), where both values are the predicted values. 

ccg r , inferred ( 0-6) = P Iw.meas (0-6), equals the aq r> ' (0-6) inferred by the radiometer measurement. 
aqf (0-6) = predicted value listed in Table 4. 
q total, meas = Calorimeter measurement = aq r ~ (0-18) + q c 

&q r , inferred ( 0-18) = q to tai, meas - q c , where q c is the predicted value. The values not in parentheses use the super-catalytic q c , while values in 
parentheses use the non catalytic q c . 
aq r ~ (0-18) = predicted value listed in Table 4. 

aq, - inferred (6-18) = ceq r ~ inferred (0-18) - aq r> 'inferred (0-6). The values not in parentheses use the super-catalytic value for aq r ' i „f ened ( 0-18), while values 
in parentheses use the non-catalytic value. 
aq r ~ (6-18) = predicted value listed in Table 4. 


Table 6. Comparison of the predicted and experimentally inferred values, defined in Table 5, but now accounting for the shock-wave 


curvature by multiplying ft, aq, (0-6) and aq,:( 6-18) by 0.85. 


t(s) 

Otyr, inferred 

(0-6) 

aq,: 

(0-6) 

Q^Jr, inferred 

(6-18) 

aq r ' 

(6-18) 

1634.0 

3.4 

4.5 

-45.4 (58.6) 

10.1 

1636.0 

12.9 

14.2 

-26.9 (90.0) 

36.3 

1637.5 

31.4 

29.8 

-1.6 (126) 

71.3 

1640.5 

81.7 

95.2 

133 (258) 

166.6 

1643.0 

151.0 

135.2 

141 (220) 

179.4 

1645.0 

63.9 

91.0 

259 (345) 

85.8 

1648.3 

14.0 

11.2 

-16.0 (640) 

10.5 


VII. Conclusions 

The radiative heating for the Fire II flight experiment was studied using a radiation-flowfield coupled analysis 
consisting of a newly developed thermochemical nonequilibrium VSL analysis and a non-Boltzmann radiation 
model. Past studies of the Fire II flight data were shown to provide a relatively wide scatter of values, none of which 
compared well with the flight data in both the nonequilibrium and equilibrium regions of the trajectory. The present 
analysis was shown to provide a good comparison with the data in both of these trajectory regions. The updated 
atomic-line and non-Boltzmann models are the reason for this good comparison relative to previous models. The 
influence of chemical equilibrium and the Boltzmann assumption of the electronic state populations were studied. 
These were shown to have only a small influence on the radiative heating at peak heating conditions and a relatively 
large influence at the earlier nonequilibrium conditions. The influence of radiation-flowfield coupling was shown to 
reduce the radiative heating by about 30%, while also slightly reducing the convective heating. An analysis of the 
calorimeter data indicated that the present model provides reasonable values for the vacuum ultraviolet radiation, 
which is shown to be a significant contributor to the total radiation. 
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